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Abstract 

We measure the branching ratio of the purely leptonic decay of the meson with unprecedented 
precision as B{D+ ^i+u) = (3.82 ± 0.32 ± 0.09) x lO^'', using 818 pb^^ of data taken on the 
^(3770) resonance with the CLEO-c detector at the CESR cohider. We use this determination 
to derive a value for the pseudoscalar decay constant /£>+, combining with measurements of the 
D+ lifetime and assuming \Vcd\ = We find /d+ = (205.8 ± 8.5 ± 2.5) MeV. The decay rate 

asymmetry p^^|^j^||^^~p^^_^||^_^^ = 0.08 ±0.08, consistent with no CP violation. We also set 90% 
confidence level upper limits on B(£)+ t+v) < 1.2 x lO"^ and B{D+ e+v) < 8.8 x lO"^. 

PACS numbers: 13.20.Fc, 12.38.Qk, 14.40.Lb 
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I. INTRODUCTION 

Purely leptonic decays of heavy mesons involve both weak and strong interactions. The 
weak part is easy to describe as the annihilation of the quark antiquark pair via the Standard 
Model W'^ boson; the Feynman diagram for D'^ — > i'^v is shown in Fig. 1. The strong 



interactions arise due to gluon exchanges between the charm quark and the light quark. 
These are parameterized in terms of the "decay constant" for the D'^ meson /^n- . The decay 
rate is given by [1] 



where Gp is the Fermi coupling constant, Md+ is the mass, me is the mass of the 
final state lepton, and Vcd is a Cabibbo-Kobayashi-Maskawa (CKM) matrix element [2], 
whose magnitude is set equal to 0.2256, the value of Vus [3]. Thus, within the context of 
the Standard Model (SM), measurement of this purely leptonic decay provides a means of 
determining and similarly measuring the purely leptonic decay of the Df meson allows 
us to determine fos- 

Meson decay constants in the B system are used to translate measurements of BB mixing 
to CKM matrix elements. Currently, it is not possible to determine /b accurately from 
leptonic B decays, so theoretical calculations of fs must be used. Since the Bg meson does 
not have fju decays, it will never be possible to determine fs^ experimentally, so again theory 
must be relied upon. If calculations disagree on D mesons, they may be questionable on B 
mesons. If, on the other hand new physics is present, it is imperative to understand how it 
affects SM-based predictions of the B decay constants. 

These decay constants can be calculated in theories of QCD. A recent calculation by 
Follana et al. [4] using an unqucnched lattice technique predicts = (207 ± 4) MeV and 
/d^ = (241 ± 3) MeV. The latter result differs by more than three standard deviations from 
the average of CLEO and Belle measurements [5] . 

Dobrescu and Kronfcld point oTit that this discrepancy can be caused by the presence 
of non-SM objects participating virtually in the decay [6]. They give three possibilities: 
(1) a new boson of charge +1 interfering with the SM annihilation, (2) a charge +2/3 
leptoquark, and (3) a charge minus 1/3 leptoquark. The charge +1 boson could either be a 
W'^ or a charged Higgs. They propose a specific two Higgs doublet model where one doublet 
gives the c, u and Icptons mass, but not the rf, s, 6, or t, and has a vacuum expectation value 
of about 2 GeV. Such a model predicts that the ratio of widths r(D+ — > t^i/)/T[D^ 
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FIG. 1: The decay diagram for t^v. 




(1) 
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is the same as the Standard Model expectation, which is in agreement with the CLEO 
measurements. 

The previous CLEO determination of = (222.6 ± 16. Tlg^) MeV is consistent with the 
FoUana et al. calculation at the one standard deviation level, but the experimental error was 
too large to provide a precision test. Here we provide a measurement based on a three times 
larger data sample and a ~15% larger efficiency based on improved analysis techniques. 

One other fully unquenched lattice calculation exists in the literature [7], although it has 
significantly larger errors than FoUana et al. [4]. Quenched calculations have also been per- 
formed [8-11], and other methods have been used [12-18]. The various theoretical predictions 
of fD+ range from 190 MeV to 350 MeV. Because of hehcity suppression, the electron mode 
— > e^v has a very small rate in the Standard Model [19]. The expected relative widths 
are 2.65 : 1 : 2.3 x 10~^ for the r+z/, /U+z/, and e+z/ final states, respectively. Unfortunately 
the mode with the largest branching fraction, t^v, has at least two neutrinos in the final 
state and is difficult to detect in decay. 

The CLEO-c detector is equipped to measure the momenta and directions of charged par- 
ticles, identify charged hadrons, detect photons, and determine their directions and energies 
with good precision. It has been described in more detail previously [20-23]. 

II. DATA SAMPLE AND SIGNAL SELECTION 

In this study we use 818 pb~^ of CLEO-c data collected from e^e" coUisions at the ■0(3770) 
resonance. This work contains our previous sample as a subset and supersedes our initial 
efforts [23]. At this energy, the events consist mostly of pure D+D-, D^D , three-flavor 
continuum, with small amounts of other final states such as ^/i!j{2S) and r+r^. 

We examine all the recorded hadronic events and retain those containing at least one 
charged D candidate in the modes listed in Table I. We use this sample to look for cases 
where we have only a single muon candidate whose four-momentum is consistent with a 
two-body D decay into a muon and a neutrino and no other charged tracks or excess neutral 
energy are present. Track selection, particle identification, tt", Ks and muon selection criteria 
are identical to those described in Reference [23] , with one important exception. The angular 
acceptance of the muon has been widened to cover 90% of the solid angle rather than 
81%. Muons deposit less than 300 MeV of energy in the calorimeter 98.8% of the time, 
while hadrons often interact and deposit significantly more energy. Thus, we define two 
cases in this paper, where case (i) refers to muon candidate tracks that deposit <300 MeV 
and case (ii) is for candidates depositing > 300 MeV, as was done previously for both our 
D'^ — > T'^f and — > /I'^i/ analyses [24, 25]. Briefly, we determine the efficiency on muons 
from e~^e~ ^'^fi~ events and compare with our Monte Carlo projection. The excellent 
agreement allows us to use the Monte Carlo efficiency for the lower energy muons observed 
in this analysis. Pion's deposit <300 MeV 55% of the time as determined from a relatively 
pure sample of — > K~7r'^ events, and their charge-conjugates. 

III. RECONSTRUCTION OF CHARGED D TAGGING MODES 

Tagging modes are fully reconstructed by first evaluating the difference in the energy, AE, 
of the decay products with the beam energy. We require the absolute value of this difference 
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to contain 98.8% of the signal events, i.e. to be within 
(rms) width of the peak value. The rms widths vary from 



i2.5 times the root mean square 
i7 MeV in the K^K~t:~ mode to 



il4 MeV in the K^ti tt tt^ mode. For the selected events we then view the reconstructed 



D beam-constrained mass defined as 



(2) 



where i runs over all the final state particles of the tag. Since the CESR beams have a 
crossing angle, we work in the center-of-mass frame. The beam-constrained mass has better 
resolution than merely calculating the invariant mass of the decay products since the beam 
has a small energy spread. Besides using D~ tags and searching for — > /x+i^, we also use 
the charge-conjugate D'^ tags and search for D~ — > jjTV; in the rest of this paper we will 
not usually mention the charge- conjugate modes explicitly, but they are always used. 

The distributions for all D~ tagging modes considered in this data sample are shown 
in Fig. 2. To determine the event numbers we first fit the mec distributions to a signal 
function plus a background shape. Then we use the signal shape to define the lower and 
upper limits in mBC, and count the number events above the background function within 
the limits. 

For the background we fit with a shape function analogous to one first used by the ARGUS 
collaboration [26] which has approximately the correct threshold behavior at large mRC- This 
function is 



/background ("^BC) = a{mBC + b) 



i 



, / mBc + h\ 

1 — exp 





ruBc + b 
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c 





(3) 



where a is the overall normalization and b, c, and d are parameters that govern the shape. 
To fix the shape parameters in each mode, we fit this function to data selected by using AE 
sidebands defined as 5a < \AE\ < 7.5(T, where a is the rms width of the AE distribution. 

For the signal we use a lineshape similar to that used for extracting photon signals from 
electromagnetic calorimeters, because of the tail towards high mass caused by initial state 
radiation [27]. The functional form is 



/signal (?7T.BC) 



A ■ exp 



_1 ( mBC-mo 
2 \ fiTiRn 



for mBc <mr,-a- a^^^ 
for mBc > rriD-Oi- cr^Bc 



Here A ^ = a 



+ 



(4) 

{l -I- erf (:^)) , "^bc is the measured mass, mu is the 
"true" (or most likely) mass o-mec mass resolution, and a and n are shape parameters. 

Table I lists the modes along with the numbers of signal events and background events 
within the signal region defined as containing 98.8% of the signal events with mec below the 
peak and 95.5% of the signal events above the peak. 

We retain the events within the mass cuts illustrated in Fig. 2 for further analysis. This 
sample includes 460,055±787±2,760 signal events, where the last error is systematic. Because 
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FIG. 2: Beam-constrained mass distributions for different fully reconstructed D~ decay candidates 
in the final states: (a) K'^-k^t:^ , (b) K^tt~tt^tt^ , (c) Kst^~ (d) iCs'7r~7r^7r^, (e) i^57r~7r'^, and (f) 
K^K^TT^ . The solid curves show the sum of signal and background functions. The dashed curves 
indicate the background fits. The region between the arrows contains the events selected for use in 
further analysis. 
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TABLE I: Tagging modes and numbers of signal and background events determined from the fits 
sliown in Fig. 2. 



Mode 




Signal Backj 


ground 


iiT+TT- 


"TT 


224, 778 ± 497 


5,957 


K+T^- 




71, 605 ± 359 


37,119 


Ks7r~ 




32, 696 ± 189 


1,576 


KsTr~ 


7r~7r"'" 


52,554 ±315 


26, 352 


Ks7r~ 




59, 298 ± 289 


14,837 


K+K 




19,124±159 


3,631 


Sum 




460, 055 ± 787 


89,472 



of their low multiplicity, it is easier to find tags in simple ji^v events than in typical D^D~ 
events. Therefore, to calculate the branching fraction we increase the number of tags by 
(1.54±0.36)%, as determined by Monte Carlo simulation. The systematic error on the signal 
number includes this uncertainty added in quadrature with the change given by varying the 
background function. 

IV. £>+ fi+v SELECTION CRITERIA 

Using our sample of D~ event candidates we search for events with a single additional 
charged track presumed to be a Then we infer the existence of the neutrino by requiring 
a measured value of the missing mass squared (MM^) near zero (the neutrino mass), where 

MM^ = (Ebeam ' E^+f - (-p^- - p^+)' , (5) 

here p/j- is the three-momentum of the fully reconstructed D~ , and £^^+(p^+) is the energy 
(momentum) of the candidate /x"*". 

To restrict the sample to candidate ^'^v events resulting from the other we exclude 
events with extra neutral energy, or more than one additional track with opposite charge to 
the tagged D, which we take to be the muon candidate. We allow such extra tracks if their 
distance of closest approach from the beam collision point is outside a region more than 5 
cm along the beam or more than 5 mm perpendicular to the beam; we do not wish to veto 
these tracks as they are usually due to interactions of the tracks from the tagging D~ in 
the calorimeter. We reject events with extra fully reconstructed Ks ti^ti~ candidates. 
We also veto events having a maximum neutral energy cluster of more than 250 MeV. This 
criterion is highly effective in reducing backgrounds especially from — > tt'^tt^ decays. We 
consider only those showers that do not match a charged track within a connected region. A 
connected region is a group of adjacent crystals with finite energy depositions. This reduces 
the probability of a false veto due to hadronic shower fragments that would otherwise show 
up as unmatched showers. 

Sometimes the decay products of the tagging D~ interact in the detector material, mostly 
the EM calorimeter, and spray tracks and neutral energy back into the rest of the detec- 
tor. We evaluate the size of these contributions to the inefficiency caused by imposing the 
250 MeV extra neutral energy requirement by using fully reconstructed D^D~ events. We 
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start with events where the K'tt'^tt'^ and the D" — > K'^tt'tt'. We then look for 

extra photons with energies >250 McV. This measures the square of the efficiency for the 
case of K~n~^7r'^ tags, our largest mode. We then measure the inefficiency for each tag mode 
by looking for fully reconstructed events where one D decays into K^tt^tt^ and the other 
into one of the other tag modes. The weighted average over all our tag modes gives an 
efficiency for our extra energy veto of (95.9±0.2±0.4)%. The details are given in Appendix 
A. 

We define 6 as the angle with respect to the positron beam direction. The muon candidate 
direction is required to have | cos^| < 0.90, and deposit less than 300 MeV of energy in the 
calorimeter, characteristic of a minimum ionizing particle. 

The MM^ from Monte Carlo simulation is shown in Fig. 3 for the proper mix of tag 
modes. The signal is fit to a sum of two Gaussian distributions with the wider Gaussian 
having about 30% of the area independent of tagging mode. The average resolution (cr) is 
defined as 

(T = fiai + (1 - fi)a2, (6) 

where ai and (72 are the individual widths of the two Gaussians and /i is the fractional 
area of the first Gaussian. The resohition of 0.0266±0.0006 GcV^ is consistent among all 
the tagging decay modes when restricting the fit range to —0.2 <MM^ <0.2 GcV^. In a 
narrower range, -0.1 kMW^ < 0.1 GeV^, the resolution is a = 0.0248 ± 0.0006 GeV^. We 
use differences in the signal function width to evaluate the systematic error. 
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FIG. 3: Monte Carlo simulation of the MM^ distributions for ^i^v events opposite the 

proper mixture of tag final states. The fit is to two Gaussian distributions centered at zero where 
the second Gaussian constitutes around 30% of the area. 

We check our simulations by using the — > Ks7r~^ decay. Here we choose events with 
the same requirements as used to search for /x+z/ but require one additional found Ks- The 
MM^ distribution for this final state is shown in Fig. 4(a) and peaks as expected at the Ks 
mass-squared of 0.25 GeV^. The resolution depends slightly on the fitting range, which must 
be specified since the data have a high MM^ background. In the interval 0.05 <MM^ < 0.35 
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FIG. 4: MM^ distribution for the decay — > Kgir^ from data (a), and signal Monte-Carlo 
simulation (b). 



GeV^ the data show a resolution oi a = 0.0247 ± 0.0012 GeV^ while the Monte Carlo fit 
gives a consistent value of ci = 0.0235 ± 0.0007 GeV^. 

The MM^ distributions for our tagged events requiring no extra charged tracks besides 
the muon candidate and no extra showers above 250 MeV as described above are shown 
in Fig. 5. We see a peak near zero mostly due to the — >■ fx^v mode we are seeking. 
The large peak centered near 0.25 GeV^, far from our signal region, results from the decay 
'K^tt^, and is expected since many Kl escape our detector. 



V. BACKGROUND CONSIDERATIONS 

In this section we will estimate backgrounds from specific sources and also specify shapes 
of several distinct background distributions. Our procedure will be to fit the signal MM^ 
distribution with the sum of the signal and background shapes and then subtract off any 
residual backgrounds, which we will show are very small. The signal shapes include both 
the ii'^u and r+i/, t+ —>■ tt+i/ distributions, separately. 

There are several background sources we need to evaluate. These include background from 
other modes, background from misidentified D^D events and continuum background 
including that from e+e~ '-filj{2S), termed "radiative return." Hadronic sources need to 
be considered because the requirement of the muon depositing less than 300 MeV in the 
calorimeter, while 98.8% efficient for muons, rejects only 45% of pions. 

We include a calculated background from — > 7r+7r° in the fit, both the shape and the 
normalization. This mode is the most difficult to reject because the MM^ peaks very close 
to zero, at 0.018 GeV^, well within our resolution of 0.0266 GeV^. It is possible for the 
photons from the tt^ decay to inadvertently be matched to the tracks from the tagging D~ 
or be missed, even though at least one photon from the 7r+7r° mode exceeds our 250 MeV 
calorimeter energy requirement and should in most cases cause such a decay to be vetoed. 
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FIG. 5: MM^ using D tags and one additional opposite sign charged track and no extra energetic 
showers (see text). The insert shows the signal region for n^v on a log scale. 



Both the shape in MM^ and the rate are accurately determined [28]. Using Monte Carlo 
simulation, we find efficiencies of 1.53% and 1.06%, for the calorimeter energy deposition 
cases (i) and (ii). respectively. (Recall case (i) is for energies less than 300 MeV, and case 
(ii) for larger energy depositions.) Multiplying this efficiency by the number of tags and 
branching ratio, (1.3=b0.2)%, gives a 9.2 event background. The uncertainty in the branching 
ratio is included in the systematic error. 

The K TT"*" mode gives a large peak in the MM^ spectrum near 0.25 GeV^. While it is 
many standard deviations from our signal region, we need to know the shape of the tail of 
this distribution. We also need to see if there are any "pathological" events due to non- 
Gaussian effects. We use the double tag events where both L>'s decay into K^tt"^ to 
evaluate both effects. Here we gather a sample of single tag K~7r'^ decays using strict AE 
and TObc criteria, and look for events with only two oppositely charged tracks where the ring 
imaging Cherenkov system (RICH) identifies one as a and other as a 7r~. The kaon is 
required to be in the RICH solid angle but the pion can be anywhere within | cos(6')| < 0.9, 
and then we ignore the kaon. The MM^ distribution is shown in Fig. 6. 

The fit gives us a rather good description of the shape of the k'^tt'^ peak, especially on 
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FIG. 6: The MM^ from events with — s- K~'ir'^ tag and the other D decaying into two tracks, 
most Hkely — > K^tt~, where the kaon is ignored. The kaon peak is fit to a double Gaussian 
distribution, containing 2,547 events. The other curve shows the expected shape for tt'^tt^ . 

the low MM^ side, where the K^tt'^tt^ background is absent. There are 2,547 Kit events. 
The small numbers of residual events peaking near the pion mass squared could be due to 
7r+7r~ events where the RICH was fooled. The fake rate in the RICH has been well measured 
as (1.2 ± 0.4)% for pions faking kaons in the momentum region of interest (see Appendix 
B). The relative branching is B{D° 7i+7i-)/B{D^ K+ti') 3.59% [2]. Thus we expect 
1.1 n~^7r~ events. There are three observed events consistent with being in the signal region 
near MM^ of zero GeV^. These three events then can be either background associated with 
K~^n~ events or tt^tt^ events. Our best estimate is that 1.9 of them are background. By 
normalizing the background based on the number of K 7r+ events in the MM^ spectrum, we 
expect 1.3 events as the background from this non-Gaussian effect in case (i) events. 

The only significant non-/i+z/ population in the signal region arises from — > r+z/. 
Out of 10,000 simulated events with D~ tags, we find events in the ^'^u signal region only 
when r+ n^u. Because of the small D^-t^ mass difference, the r"*" is almost at rest 
in the laboratory frame and thus the vr"^ has relatively large momentum causing the MM^ 
distribution to populate preferentially the low MM^ region, even though there are two missing 
neutrinos in this case. Thus, we generate a shape from Monte Carlo specifically for this one 
decay sequence as shown in Fig. 7. 

Other backgrounds from r"*" decays include additional missing particles. We form a shape 
consisting of a sum of the r"*" decay modes p^z/, /i+z/z/ and other "similar" modes p+7r°, and 
TT^lx^v. All the relevant branching ratios are known, where we take the r^v rate by 

multiplying our previous n^u result by 2.65, the Standard Model prediction. We use this 
shape to describe these backgrounds; we do not, however, fix the normalization in the fit. 
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FIG. 7: The simulated MM^ from — > t'^u, t+ — > tt'^u fit to the sum of two Gaussian distribu- 
tions, whose widths are allowed to vary on both sides of their respective maximums. 

We have also checked the possibility of other D^D~ decay modes producing background 
with an equivalent 1.7 Po^^ Monte Carlo sample; we find no additional events. The D^D and 
continuum backgrounds are also evaluated by analyzing Monte Carlo samples corresponding 
to 4.1 and 3.0 fb~^, respectively. To normahze our Monte Carlo events to our data sample we 
used CTj^Q-^ — 3.7 nb [29] and (Tcontinuum = 18 nb. We also found no events in our analysis of 
a simulated radiative return sample equivalent to 2.7 fb~^. Our total additional background 
is 2.4±1.0 events, with the individual components listed in Table II. 

TABLE II: Backgrounds from additional sources, not contained in the fitting functions. 



Mode 


# of Events 


Continuum 


0.8±0.4 


-T^ + 

K TT^ 


1.3±0.9 


modes 


0.3±0.3 


Sum 


2.4±1.0 



VI. BRANCHING RATIO AND DECAY CONSTANT 

We preform a binned maximum liklihood fit to the case (i) MM^ distribution up to a 
MM^ of 0.28 GeV^. Beyond that value other final states such as r;7r+ and K^tt^t:^ begin to 
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contribute. The fit shown in Fig. 8 contains separate shapes for signal, tt+tt , K tt"*", t'^v 
{t^ — > 7r''"P), and the background cocktail described above. Here we assume the Standard 
Model ratio of 2.65 for the ratio of the T^vj^i^v component and constrain the area ratio 
of these components to the product of 2.65 with i3(r+ 7r+z/) = (10.90±0.07)% [2] and the 
55% probabiUty that the pion deposits <300 MeV in the calorimeter. The normalization of 
the tt+tt'^ component is also fixed at 9.2 events, the product of the number of tags, times the 
branching fraction, times the 1.53% detection efficiency. The normalization of the additional 
background shape described above is allowed to float. 

The fit yields 149.7±12.0 \x^v signal events and 25.8 r^z/, r"*" — >■ tt+z/ events (for the entire 
MM^ range). We can also perform the fit allowing the r^v^ tt'^P component to float. 

(See Fig. 9.) Then we find 153.9±13.5 fi'^u events and 13.5±15.3 r+z/, r"*" n^u events, 
compared with the 25.8 we expect in the Standard Model. Performing the fit in this manner 
gives a result that is independent of the SM expectation of the — > t^v rate. To extract 
a branching fraction, in either case, we subtract off the 2.4±1.0 events determined above to 
be additional backgrounds, not taken into account by the fit, and divide by the product of 
the efficiency and the number of tags. 

The detection efficiency of 81.8% for the single muon includes the tracking and particle 
identification efficiencies, the probability of the crystal energy being less than 300 MeV, and 
the 95.9% efficiency of not having another unmatched shower in the event with energy greater 
than 250 MeV; the latter is determined from the data presented in Table V of Appendix A. 
The systematic errors on the branching ratio are listed in Table III. 

TABLE III: Systematic errors on the _D+ fi'^u branching ratio. 



Systematic errors (%) 



Track finding 


0.7 


PID cut 


1.0 


MM^ width 


0.2 


Minimum ionization cut 


1.0 


Number of tags 


0.6 


Extra showers cut 


0.4 


Radiative corrections 


1.0 


Background 


0.7 


Total 


2.2 



The systematic error on the MM^ fit is determined by changing the signal shape and 
the fitting range. The difference in signal shapes between the Ksn^ data and Monte Carlo 
is 0.0012±0.0014 GeV^. We refit the case (i) data while increasing a by 0.0012 GeV^ 
0.0024 GcV^, and finally letting a float. (We flx the r^v/ ji^v ratio.) The resulting numbers 
of events change from our baseline by +0.41, +0.79 and +0.26 events, respectively. This 
allows us to set the 0.2% systematic error from this source. 

The track flnding and particle identiflcation efficiencies associated with the single muon 
are determined by comparing selected samples formed using partial reconstruction [30] to the 
Monte Carlo simulation. We include the particle identification because we do veto identified 
kaons as muon candidates. 



13 



1630608-012 
\ 1 




MM^ (GeV^) 

FIG. 8: Fit to the MM^ for case (i). Here the ratio of t^i^, t+ tt^u to fi^i^ events is fixed to 
the SM value. The points with error bars show the data. The black (dashed) curve centered at 
zero shows the signal fj.^u events. The dot-dashed (red) curve that peaks around 0.05 GeV^ shows 
the — > r+i/, r"*" vr+P component. The solid (blue) Gaussian shaped curve centered on the 
pion-mass squared shows the residual vr'^vr'^ component. The dashed (purple) curve that falls to 
zero around 0.03 GeV^ is the sum of all the other background components, except the K vr"*" tail 
which is shown by the long-dashed (green) curve that peaks up at 0.25 GeV^. The solid (black) 
curve is the sum of all the components. 
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FIG. 9: Fit to the MM^ for case (i) allowing the t^u, t+ — > ir^u component to float. The points 
with error bars show the data. The black (dashed) curve centered at zero shows the signal ^u^i^ 
events. The dot-dashed (red) curve that peaks around 0.05 GeV^ shows the — r+z^, t+ — > vr+P 
component. The solid (blue) Gaussian shaped curve centered on the pion-mass squared shows the 
residual vr'^vr'^ component. The dashed (purple) curve that falls to zero around 0.03 GeV^ is the 
sum of all the other background components, except the K tt^ tail which is shown by the long- 
dashed (green) curve that peaks up at 0.25 GeV^. The solid (black) curve is the sum of all the 
components. 
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MM2 (GeV^) 

FIG. 10: Fit to the MM^ for case (ii) which has httle fi'^u signal contribution, and tests our 
understanding of the background. The points with error bars show the data. The fi^u component, 
the black (dashed) curve (almost invisibly small) shows the signal /x"*"!^ events fixed from the case 
(i) fit. The points with error bars show the data. The dot-dashed (red) curve that peaks around 
0.05 GeV^ shows the t^v, ir^u component. The solid (blue) Gaussian shaped curve 

centered on the pion-mass squared shows the residual vr^'^vr'^ component. The dashed (purple) curve 
that falls to zero around 0.03 GeV^ is the sum of all the other background components, except the 
K Tr"*" tail which is shown by the long-dashed (green) curve that peaks up at 0.25 GeV^. The solid 
(black) curve is the sum of all the components. 
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A check of the background is provided by considering case (ii) , where more than 300 MeV 
is deposited in the calorimeter by the muon candidate track. Only 1.2% of muons pass such 
a requirement. We fit this sample as in case (i), but here fixing both the fi^u and the r+z/ 
contributions from the case (i) fit (with the ratio of the two fixed). The normalizations of 
the K^n'^ tail and the background shape are allowed to float. The fit is shown in Fig. 10. 
The number of events in the signal region, MM^ <0.05 GeV^, is 1.7, fixed from the /j.'^i' 
sample, 5.4 fixed from the 7r+7r°, and 4.0 from the r+z/. This sums to 11.1 events, while we 
count 11 events in this region. Thus we have an excess of -0.1±3.3 events, which is consistent 
with our other background estimate of 2.4±1.0 events and gives us confidence in using this 
estimate. 

The branching fraction determined from fixing the r'^u contribution relative to the /J^'^i', 

is 

B{D+ fi+u) = (3.82 ± 0.32 ± 0.09) x 10"^ . (7) 

The decay constant f£,+ is then obtained from Eq. (1) using 1040±7 fs as the lifetime 
[2] and 0.2256 as l^cdl [3]. Our final result is 

= (205.8 ± 8.5 ± 2.5) MeV . (8) 

A somewhat less precise value is obtained by fioating the r+i^ to fi'^v ratio. That fit gives 

B{D+ ^+i/) = (3.93 ± 0.35 ± 0.09) x 10"^ . (9) 

The corresponding value of the decay constant is 

fjj+ = (207.6 ± 9.3 ± 2.5) MeV . (10) 

The former value is the most precise measurement in the context of the Standard Model, 
while the latter does not use any Standard Model assumptions. In both cases the additional 
systematic errors due to the lifetime measurement and the error on \Vcd\ = \Vus\ are 
negligible. 

The data have already been corrected for final state radiation of the muon, as our Monte 
Carlo simulation incorporates this effect [31]. There is however, another process where 
the '-jD*^ 7Ai^^, where the D*^ is a virtual vector or axial-vector meson. The 

D*^ — > ii'^v transition is not helicity-suppressed, so the factor a for radiation is compensated 
by a relative factor {Mj)+ /m^Y. Using Eq. (12) of Burdman et al. [32] and imposing the 
250 MeV photon cut, we find that the radiative rate is approximately 1%, to which we 
assign a ±1% systematic error. This is essentially the same calculation done by Dobrescu 
and Kronfeld for Df — > /i+z/ decays [6]. (The results shown above for the branching fractions 
and /d+ are all radiatively corrected; the branching fractions have been reduced by 1%.) 

YII. SEARCH FOR D+ t+u 

We also use our data to perform a search for the r'^u final state. Here we do a simultaneous 
binned maximum liklihood fit to both the case (i) and case (ii) data fixing the ratio of the 
T^u final state to be 55/45 in the two cases, determined by the relative acceptances for the 
300 MeV calorimeter energy requirement. The fits are shown in Fig. 11. 
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FIG. 11: Fit to the MM^ for cases (i) (top) and (ii) (bottom) with the t^v components fixed in 
the ratio 55/45. The points with error bars show the data. The black (dashed) curve centered at 
zero shows the signal fj.~^i> events. The dot-dashed (red) curve that peaks around 0.05 GeV^ shows 
the — > r"*"zv, r"^ 7r"'"P component. The solid (blue) Gaussian shaped curve centered on the 
pion-mass squared shows the residual vr'^vr'^ component. The dashed (purple) curve that falls to 
zero around 0.03 GeV^ is the sum of all the other background components, except the K tt'^ tail 
which is shown by the long-dashed (green) curve that peaks up at 0.25 GeV^. The solid (black) 
curve is the sum of all the components. 



The fit yields a sum of 27.8±16.4 t^u, r+ — > tt^u events for the entire MM^ range. To be 
conservative in setting an upper limit, we assume all events are signal and do not subtract 
additional backgrounds from this yield. We include the small systematic errors from the 
fitting procedure in our calculations. We find 

}3{D+ T^v) < 1.2 X 10"^ (11) 

at 90% confidence level, and the ratio to the rate divided by the Standard Model 
expectation of 2.65 is 

2.65 • T{D+ ^ /i+z/) ^ ' 
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also at 90% confidence level. 



VIII. SEARCH FOR D+ DECAY INTO AN POSITRON PLUS NEUTRINO 

We use the same tag sample. Candidate positrons are selected on the basis of a likelihood 
ratio constructed from three inputs: the ratio between the energy deposited in the calorimeter 
and the momentum measured in the tracking system, the specific ionization dE/dx measured 
in the drift chamber, and RICH information. Other criteria remain the same, except that 
wc require that the positron candidate track be in the calorimeter barrel with | cos^| < 0.81. 
We do not find any candidates allowing us to set a limit 

^ e+i/) < 8.8 X 10-^ at 90% c.L, (13) 
which is three orders of magnitude above the SM prediction. 

IX. CONCLUSIONS 

The result shown here represents the only precision measurement of the pscudoscalar 
decay constant //)+. We have significantly improved our previous results. The statistical 
error has been reduced by almost a factor of two. The systematic errors remain small. This 
result uses all the CLEO-c data collected at the ■?/'(3770) and, as such, all previous results 
are superseded. 

The branching fraction, assuming the Standard Model ratio for T^v/jji^v is 

B{D+ /i+z/) = (3.82 ± 0.32 ± 0.09) x 10"^ , (14) 
and the decay constant is 

/^+ = (205.8 ± 8.5 ± 2.5) MeV . (15) 
If, on the other hand, we allow the t^u contribution to fioat, we find 

B{D+ li+v) = (3.93 ± 0.35 ± 0.09) x 10"^ . (16) 
The corresponding value of the decay constant is 

/^+ = (207.6 ± 9.3 ± 2.5) MeV . (17) 

These results are all radiatively corrected. 

Our new values are consistent with our previous measurement [23] , as well as the upper 
limit set by Mark III [33], and the results based on reported yields of 1 and 2.7 events from 
BES I and II [34], respectively. We also determine fos/ fD+ = 1-326 ±0.075, using the world 
average value of absolute measurements for — > £+z/ as compiled by Rosner and Stone [5], 
where we include the radiative correction also on the Df rate. 

Our result for /£>+ , is consistent with the most accurate unquenched lattice QCD (LQCD) 
calculation of Follana et al. who give a value of (207±4) MeV. This implies that the somewhat 
greater than three standard deviation discrepancy of the experimental measurements of 
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/^^ = (273 ± 10) MeV [5] with the FoUana et al. prediction of (241±3) MeV cannot be 
explained by how they handle the charm quark in their calculation. In fact, since the s 
quark is heavier than the d quark, it should be easier for lattice calculations to predict 
than [35]. It may be the case that physics beyond the Standard Model is raising the value 
of fo^ in one of the ways suggested by Dobrescu and Kronfeld [6], or via R-parity violating 
supersymmetry [36] . Other Standard Model based predictions are listed in Table IV. 



TABLE IV: Theoretical predictions of /£>+ and /^+ / fu+ 



Theory 


/d+ (MeV) 


fD+/fD+ 


LQCD (HPQCD+UKQCD) [4] 


207 ±4 


1.164 ±0.011 


LQCD (Fermilab+MILC) [7] 


201 ± 3 ± 17 


1.24 ±0.02 ±0.07 


QL (QCDSF) [8] 


206 ± 6 ± 3 ± 22 1.07 ± 0.02 ± 0.02 


QL (Taiwan) [9] 


235 ± 8 ± 14 


1.13 ±0.03 ±0.05 


QL (UKQCD) [10] 


210 ± 10t\l 


1.13±0.02l|];[]^ 


QL [11] 


211 ± 141^2 


1.10 ±0.02 


QCD Sum Rules [12] 


177 ±21 


1.16 ±0.01 ±0.03 


QCD Sum Rules [13] 


203 ± 20 


1.15 ±0.04 


Field Correlators [14] 


210 ± 10 


1.24 ± 0.03 


QCD Sum Rules [15] 


195 ± 20 




Relativistic Quark Model [16] 


234 


1.15 


Potential Model [17] 


238 


1.01 


Isospin Mass Splittings [18] 


262 ± 29 





It is possible in some models of new physics that there is a difference in the jiu decay rate 
between and D~ mesons, due to a CP violating interaction [37]. Separating our data 
into these two classes we find 228,945±551 tags and 231,107±552 D~ tags. Fitting the 
data by fixing the relative r^i/, — > n^u contribution relative to /x^z^, we find 76.0±8.6 
li'^u events and 64.8±8.1 /I'v events. The resulting CP violating asymmetry is 

. '2] = 0.08 ± 0.08 . (18) 

At 90% confidence level the limits are -0.05 < Acy> < 0.21. 

We do not find positive evidence of the decay L>+ t^v. Our limit is 

B{D+ T+v) < 1.2 X 10-^ (19) 

at 90% confidence level, and the ratio to the ii^v rate, divided by the Standard Model 
expectation of 2.65 is 

2.65. r(c.^^) < P«> 

also at 90% confidence level. 

Some non-standard models predict significant rates for the helicity suppressed decay 
L>+ e^v [38]. Our upper limit of 8.8 x 10~^ at 90% c.l. restricts these models. 
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APPENDIX A: DETERMINATION OF THE EFFICIENCY OF THE 250 MEV 
CRITERIA ON ADDITIONAL PHOTONS 



Although we do not expect more than a few percent inefficiency due to rejecting events 
with an additional neutral energy cluster > 250 MeV, we do not want to incur a large 
systematic error due to this potential source. Therefore we perform a full five-constraint 
kinematic fit to the double tag event samples, where one D decays into K^tx^tt^ and the 
other into one of the other tag modes. The constraints are that the total energy sum to twice 
the beam energy, the total three momentum be zero, and the invariant masses of the two D 
candidates be equal. We do not require them to equal the known mass. The result of 
this fit is a common D candidate mass and a x^. Restricting our samples to low virtually 
eliminates all backgrounds at the expense of some signal. Specifically, we require that the 
probability of x^, for five constraints be greater than 1%, which eliminates 32% of all event 
candidates. The numbers of events in the decay modes we use are hsted in Table V. 



TABLE V: Numbers of D'^D events and the efficiency for the first mode when an extra photon 
> 250 MeV is also required. 



Mode 1 


Mode 2 


Events iViost 




Mev) £250 (%) of Mode 1 


K~^TT~TT~ 




4389 


431 


95.0±0.2 






2590 


208 


96.8±0.6 






1255 


112 


95.9±0.8 






1885 


153 


96.8±0.7 


KsTT-TT^ 




2648 


205 


97.7±0.5 


K+K'TT- 




714 


75 


94.2±1.1 



Weighted Average 95.9±0.2 



To first order the fully reconstructed D'^D" — > (ir+7r~7r~) {K'n'^n'^) can be considered 
the superposition of two single tag D'^ fi'^u candidate events where the single tag is 
K~^'K~7r~. Then the efficiency of the 250 MeV cut in the Km: mode is given by 

£250'' = - Niost/NKTrn-KTzn)- (Al) 

We then combine the large and precise K'n'^Tr'^ mode with each of the other tags in turn, 
where 

^250'^'^ = (1 ~ Niost/NK-K'K-raodc)/^2^'' ■ (A2) 

This method ensures that the number of interactions of particles with material is the same 
as in the tag sample used for the ^I'^v analysis. 

The results are hsted in Table V. The numbers of events hsted are those with a x^ cut 
applied. The overall efficiency for accepting the double tag event requiring that there not 
be any photons above 250 MeV is given along with the derived efficiency for each mode. 
The weighted average over all of our tag modes is (95.9±0.2±0.4)%. The systematic error 
arises only from the consideration that we have analyzed a situation corresponding to two 
overlapping tags rather than one tag plus a muon. 
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APPENDIX B: RICH PARTICLE IDENTIFICATION EFFICIENCIES 

For two-body decays of D mesons, most of the particle identification ability in CLEO 
comes from the RICH detector. Information on the angle of detected Cherenkov photons is 
translated into a likelihood of a given photon being due to a particular particle. Contributions 
from all photons associated with a particular track are then summed to form an overall 
likelihood denoted as Ci for each particle hypothesis. To differentiate between pion and 
kaon candidates, wc use the difference: — 21og(£^) + 21og(£i^). A value of zero is used 
to distinguish between the two possibilities. We require a minimum of three Cherenkov 
photons. 

Here wc use a selected sample of D^D — > K^n^ + i^^vr^ decays. We use —0.0194 < 
AE < 0.0175 GeV and 1.8617 < mec < 1-8673 GeV for both candidates. This is essentially 
a background free sample. We expect only K~Tr~^K~^Tr~ decays since doubly Cabibbo sup- 
pressed decays are forbidden due to quantum correlations and the mixing rate as measured 
is small enough not to allow us to see any events. The momentum distribution of the tracks 
is flat between 700 MeV/c and 1 GeV/c. 

The results are shown in Table VI. The first column labeled "No ID" gives the number of 
K~n~^; K~^n~ pairs called right sign (RS) and the number of i^^7r+; K~tt~^ (or K^tt~; K^ti^) 
pairs that are wrong sign (WS) using only the kinematical constraints of /S.E and mec given 
above. The subsequent columns show the results of applying the RICH particle identification 
criterion to identify only the kaons, only the pions and then both kaons and pions. 

TABLE VI: Results of RICH identification on double tag events. RS indicates right sign and WS 
indicates wrong sign events. 

Mode No ID Single K Single vr Double ID 

(RS) (WS) (RS) (WS) (RS) (WS) (RS) (WS) 

K-TT+\ K^TT^ 1896 914 1717 11 1846 24 1673 1 



The relevant results are summarized as: 

• The pion efficiency is (97.3±0.3)%. 

• The kaon efficiency is (90.6±0.7)%. 

• The rate of pions faking kaons is (1.2 ± 0.4lo.i)%- 

• The rate of kaons faking pions is (2.6 ± 0.5io.i)%- 

The one doubly identified wrong sign event could be a mixed event, although that is 
rather unlikely. We use it to assign a negative systematic error on the fake rates in case 
there is background in our sample. 
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